SHEP 95/42 
KANAZAWA 96-02 
HUB-EP-96/6 
HLRZ 77/95 
WUB 95-45 

Dual Superconductor Scenario of Confinement: 
A Systematic Study of Gribov Copy Effects 



G.S. Bali 

Physics Department, The University, Highfield, Southampton SO 17 IB J, UK 

V. Bornyakov 

Department of Physics, Kanazawa University, Kanazawa 920-11, Japan 
and Institute for High Energy Physics, 142284 Protvino, Russid^\ 

M. Miiller-Preussker 

Institut fur Physik, Humboldt Universitdt zu Berlin, 10099 Berlin, Germany 

K. Schilling 

Fachbereich Physik, Bergische Universitdt, Gesamthochschule, 42097 Wuppertal, Germany 
and HLRZ c/o KFA, 52425 Julich, and DESY, 22603 Hamburg, Germany 

Abstract 

We perform a study of the effects from maximal abelian gauge Gribov copies in the 
context of the dual superconductor scenario of confinement, on the basis of a novel 
approach for estimation of systematic uncertainties from incomplete gauge fixing. 
We present numerical results, in SU (2) lattice gauge theory, using the overrelaxed 
simulated annealing gauge fixing algorithm. We find abelian and non-abelian string 
tensions to differ significantly, their ratio being 0.92(4) at (3 = 2.5115. An approxi- 
mate factorization of the abelian potential into monopole and photon contributions 
has been confirmed, the former giving rise to the abelian string tension. 
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1 Introduction 



Understanding how confinement arises from QCD is a central problem of strong interac- 
tion physics. G. 't Hooft and S. Mandelstam proposed the QCD vacuum state to behave 
like a magnetic superconductor |], 0. A dual Meissner effect is believed to be responsible 
for the formation of thin string-like chromo-electric flux tubes between quarks in SU(N) 
Yang-Mills theories. This confinement mechanism has indeed been established in compact 
QED || [|, ||. The disorder of the related topological objects — magnetic monopoles - 
gives rise to an area law for large Wilson loops and, thus, leads to a confining poten- 
tial. The application of this idea to non-abelian gauge theories is based on the abelian 
projection ||, reducing the non-abelian SU(N) gauge symmetry to the maximal abelian 
(Cartan) subgroup {/(l)-^ -1 by gauge fixing the off-diagonal components of the gauge field. 
Then the theory can be regarded as an abelian gauge theory with magnetic monopoles 
and charged matter fields (quarks and off-diagonal gluons). The dual superconductor 
idea is realized if these abelian monopoles condense. In this scenario large distance (low 
momentum) properties of QCD are carried by the abelian degrees of freedom (abelian 
dominance). 

Nonperturbative investigations of this conjecture became possible after formulating 
the abelian projection for lattice gauge theories 0. In previous numerical studies (see 
reviews [R Q and references therein), it has been demonstrated that the maximal abelian 
(MA) gauge was a very suitable candidate for lattice investigations of the 't Hooft- 
Mandelstam confinement scenario, being the only known renormalizable abelian gauge. 

These studies provided strong evidence for the QCD vacuum indeed to act like a dual 
superconductor. In Ref. |L0| the value of the abelian string tension, i.e. the slope of the 
static potential between charge one static sources at large distance, has been observed to 
be close to that of the potential between static quarks in the fundamental representation 
of the non-abelian theory. This feature supports abelian dominance as predicted by 't 
Hooft. Recently, various results in favor of this picture have been obtained by other 
groups as well ||. 

At this stage it appears to be important to address such issues on a more quantitative 
level. This requires a careful study of the problem of gauge (or Gribov) ambiguities |TT| ], 
and the resulting biases on abelian observables. In the present paper we aim at remov- 
ing this uncertainty of previous studies. We shall develop a new effective gauge fixing 
algorithm, thus reducing gauge fixing ambiguities with respect to the standard relaxation 
algorithm employed so far. We propose a numerical procedure to estimate the remaining 
gauge fixing biases. This enables us to carry out high precision measurements of the 
abelian string tension and other abelian observables, with control over systematic errors. 
The main results have been obtained on a 32 4 lattice at /3 = 2.5115. 
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Our investigation revealed that, compared to typical statistical errors, the effect of 
gauge copies cannot be neglected, even with the improved gauge fixing algorithm. Never- 
theless, the new algorithm does reduce the variance of observables with respect to various 
gauge copies considerably, compared with the traditional relaxation plus overrelaxation 
algorithm, and, in average, yields larger values of the functional to be maximized. 

Our paper is organized as follows. In Section 2 we shortly review the status of gauge 
fixing and monopole kinematics to set the stage and to present the underlying formulae in a 
self contained fashion. In Section 3, we present the overrelaxed simulated annealing gauge 
fixing algorithm and compare it to the standard algorithm. In Section 4 we elaborate 
on a procedure to estimate the effect of remaining gauge fixing ambiguities on abelian 
observables, such as the monopole density or the abelian potential. Physical results on 
the non-abelian and abelian potentials, the decomposition of the potential into monopole 
and photon contributions, the abelian potential between charge two static sources and 
the monopole density are presented in Section 5. 



2 Physics from maximal abelian gauge 
2.1 Abelian projection 

We start from the lattice version |7j of 't Hooft's abelian projection ||. The idea is to 
partially fix gauge degrees of freedom such that the maximal abelian (Cartan) subgroup 
(f/ 7V ~ 1 (l) for SU(N) gauge groups) remains unbroken. 

A few abelian gauges have been suggested in Ref. |J where MA gauge, referring to 
a differential gauge condition, has been favored. Lattice simulations have indeed demon- 
strated MA gauge to be very suitable for investigations of abelian projections of gauge 
theories. 

In the following, we will restrict ourselves to the case of SU(2) gauge theory. Fixing 
MA gauge on the lattice amounts to maximizing the functional (V = N sites ), 

^) = ^E Tr K/3^) , (i) 

" * n,fj, 

with respect to local gauge transformations, 

Condition Eq. ([!]) fixes (besides other possible degeneracies) g n only up to multiplications 
g n -» v n g n with v n = exp(ia n T 3 ), r 3 = er 3 /2, -2tt < a n < 2ir, i.e. g n G SU(2)/U(1). 
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It has been shown that the corresponding continuum gauge, defined by the bilinear 
differential gauge condition, 

(dp±igA$At = , A$ = Al±iAl , (3) 

is renormalizable |12|], a feature, that is crucial for a continuum interpretation of lattice 
results. 

After a configuration has been transformed to satisfy the MA gauge condition, a coset 
decomposition is performed, 

Un,fi C|i,/iKi,(j j (4) 

where V n ^ = exp(i0 niM r 3 ), — 2-n < n)/U < 2n, transforms like a (neutral) gauge field 
and C„ jA1 like a charged matter field with respect to transformations within the residual 
abelian subgroup, 

Quark fields are also charged with respect to such £7(1) transformations. The abelian 
lattice gauge fields V n ^ constitute an abelian projected configuration. 

The SU{2) action of the original gauge theory can be decomposed into a £7(1) pure 
gauge action, a term describing interactions of the £7(1) gauge fields with charged fields, 
i.e. the off-diagonal components, and a self-interaction term of those charged fields |13| . 
Maximizing the diagonal components of all gauge fields with respect to the off-diagonal 
components amounts to enhancing the effect of the pure £7(1) gauge part in comparison 
with those contributions containing interactions with charged fields. On a heuristic level, 
one might expect the importance of the £7(1) degrees of freedom to be enhanced in the 
MA projected theory, in comparison with other abelian projections. 

The abelian Wilson loop for charge one static sources is defined as 

W a \C) = l-Tr ( I] VA = Re ( J] tti ) , (6) 

A \iec J \iec J 

where u ntfJi = exp(i8 ntlM ) , 9 n ^ = ^4> n ,^- C denotes a closed contour. 

In what follows, we will use 9 n ^ to specify abelian lattice fields for the sake of convenience. 
The "Re" symbol can be omitted from the definition of the abelian Wilson loop since 
expectation values of this operator become automatically real due to charge invariance 
(Im(n, eC ^) = 0). 

2.2 Monopoles on the lattice 

One defines magnetic monopoles with respect to the residual £7(1) gauge group in the 
way proposed in Ref. [frjj for £7(1) lattice gauge theory: abelian plaquette variables, 
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9n,[iu = 9n,iJ, + On+[i,v — 9n+u,p — 9n,v , @n,fJ,v ^ [ 47T, 471") , (7) 

can be decomposed into a periodic (regular) part, — 7r < niAll/ < 7r, and a singular part, 
m n> ^ = 0,±1,±2, 

@n,p,v @n,p.v 27lTTl n ^ v . (8) 

Qn,tiv describes the [/(l)-invariant "electromagnetic" flux through the plaquette and m n ^ u 
is the number of Dirac strings passing through it. Magnetic monopole currents k n fJj , 
residing on the links of the dual lattice, are defined as 

kn,u ~. n,vp<j9u@n,pa j (9) 
47T 

where a lattice forward derivative is used, d„f n = f n+ c — f n . It is obvious, that k n ^ 
represents a conserved current: 

E^n,M = " ( 10 ) 
A* 

Note, that the shifted currents j njA1 = k n+ ^ form closed loops on the dual lattice. As 
it has already been mentioned, first lattice results on abelian dominance have been ob- 
tained in Ref. where the string tension computed from abelian Wilson loops after MA 
projection, Eq. (|B|), has been numerically found to agree approximately with the string 
tension as extracted from the full, non-abelian theory. Since then additional lattice results 
in support of the 't Hooft-Mandelstam confinement scenario have been found || ||. In 
particular the monopole contribution to the abelian string tension has been investigated. 



2.3 Photon and monopole dynamics 

It is known that in the Villain formulation of compact U(l) lattice gauge theory one meets 
an exact factorization of expectation values of Wilson loops into monopole and photon 
parts and one would expect to encounter remnants of this in other lattice formulations 
of the theory, in the form 

(W u{1) ) w (W mon ) (W ph ) . (11) 

This factorization would induce a decomposition of the potential between static charges 
at separation R, 

V U{1) (R) w V mon (R) + V ph (R) , (12) 

where the photon contribution is expected to be Coulomb-like. 

The decomposition rule, Eq. (|l"2"D, has indeed been verified to hold approximately 
in simulations [T5j for the case of compact U(l) with Wilson action. Similarly, in MA 
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projected SU{2) gauge theory, abelian Wilson loop potentials were found to decompose 
qualitatively |16|, |17| . 

As we aim at more quantitative conclusions about the role of monopoles and photons 
within gluodynamics we will proceed to collect all necessary formulae, following essentially 
Kef. 1- 

The abelian Wilson loop operator Eq. @ can be easily transformed to 

W a \C) = exp J] 6 n ^M n ^) , (13) 

where M n ^ v is an integer valued antisymmetric field, living on plaquettes, satisfying the 
condition, d~M n ^ v = J n y, J njV being an external current associated to the (oriented) 
contour C where it takes the values ±1. d~ denotes the lattice backward derivative, 
dpfn — fn~ fn-p,- Let us now rewrite 9 n ^ in terms of dual potentials p n M , and a photon 
field 0^, 

On,iw = ^^vaf3daPn,P + dfl n ,v ~ ^n./u + ^ , (14) 

where if are zero modes defined by 



' n n 

The dual vector potential satisfies the equation, 

dud~p n ,fi - d~d u p n;V = -2nk njlM . (16) 
By imposing the Lorentz gauge condition (d^p n ^ = 0) one finds, 

p n ^ = 2n^2 D(n - m)k m ^ + const. , (17) 

m 

where D(n — m) denotes the lattice Coulomb propagator. Eqs. (|i~4"|), (|I~5|) and ([H]) define 
6' n up to an irrelevant constant. By inserting Eqs. ([14]) and (p~7|) into Eq. (|13|) , one arrives 
at 

W ab = W mon W ph W fv , (18) 

where 



W mon = exp ( -2m £ K n ^D(n - m)k m A , K Uifl = \e ViUxfi d v M n 

\ n,m / ^ 

W ph = exp (-i £ dX,„M nifl „) = exp £ ff^jj^ , 



a/3 , (19) 

(20) 



W^ = expf-^£M n ,J , (21) 

where the latter expression denotes a finite volume contribution. After applying Eq. fll^D, 
W ph can be written in a form, that is more convenient for calculations, 

W ph = exp ( -i J2 d-dn^Din - m)J m \ . (22) 

\ n,m / 



It is worth mentioning that in Refs. |L5 , 16, [17| only monopole and photon contri- 



butions have been considered while the finite volume contribution has been neglected. 
As pointed out above, for U(l) gauge theory with Villain action, it can be shown that 
Eq. ( |T8D holds for expectation values of Wilson loops as well. 

It is straightforward to generalize Eqs. (|T9|) - ( ^1]) to the case of extended mono- 
poles [P| of size I. These are defined as 



^n]ti — a ^ iivpcrd[} 9 n pa — k m ^ , (23) 

^ mGc ; 

where d®f n = f n+u> - f n , 9®^ = Emes; Qn,^, Ci and s; are cubes and squares of linear 
extent I. Subsequently, the operator W ab {C) can be decomposed in the same way as in 
Eq. (|TBD with the monopole contribution defined by 



W^ mon = e^ V {-2^K^\n-m)k^\ , K% = \e vm ,df M n ^ , (24) 

\ n,m / 

D^ l \n — m) being the lattice Coulomb propagator on a lattice of linear dimension L/l. 

It is useful to apply an alternative decomposition of the abelian Wilson loop into 
monopole and photon parts, based on the monopole vector potential ||, 6™™ = 9 n ^ — 9' n . 
From Eqs. (|8|) and ( |I~4"D one finds 9™°J l to satisfy the condition 

H'C - d " d Pn™ = 27rd;m n , pl , , (25) 

which has the Lorentz gauge solution, 

9™™ = -27rJ2D(n-m)d;m m> ^ + const. . (26) 

m 

We can now write 



W a \C) = exp ( % ]T 9T on Ji exp \i £ 9[Ji ] = W monJv (C) W ph (C) . (27) 
lec J \ lec 
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In Eq. fl2"?D, W mon 'f v combines both, monopole and finite volume contributions. This 
representation has the advantage that — provided the vector potentials 6*™° n and 9' 
have been computed — smearing techniques [jnj can be applied, which have been proven 
to be very useful in extracting potentials from Wilson loops. 



2.4 Charge two potential 



From Eq. 



the following decomposition of link variables within the adjoint represen- 



tation U„„ can be derived: 



where 



U 



C A yA 
n,n n,fi 



and 



1 



2 



la/3 



-Tr 



(TaVn^crpV^ = exp (iT 3 cj) 



T 3 is the generator of the adjoint representation and the link angle 0^ 
above. 

The "adjoint" abelian Wilson loop, 



(28) 
(29) 

(30) 

has been defined 



W ab ' adj (C) = -Ti A 
3 

can be expressed in terms of abelian link angles 



lec 



W ab,adj( C ) = (1 + 2 • W ab > 2 (C)) 

3 



with 



W ab '\C) = cos(2ft 



c 







c 



E j 



(31) 



(32) 



(33) 



The adjoint static "quark" field has one neutral and two charged components with respect 
to residual abelian gauge transformations. The neutral component does not interact with 
abelian gauge fields and gives rise to the constant on the right hand side of Eq. (|32|) . Thus, 
it is evident that, for a full description of the interaction between adjoint sources, effects 
from the off-diagonal gluon fields have to be considered. However, one still might 
expect abelian dominance to hold for charged quarks. We define the abelian projected 
analogue of the adjoint potential, V o6,a *'(i2), from the asymptotic decay of large Wilson 
loops W ab ' 2 (R,T): 



W ab ' 2 (R,T) oc exp (-V abMj (R)T) + ■ 



(34) 
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We expect the string tension extracted from V ab ' ad ^{R) to approximate the adjoint 
string tension which can be defined at intermediate distances. This expectation is sup- 
ported by an investigation of the abelian projection for the adjoint representation in the 
limit of large N c [EJO, M 



3 Efficient gauge fixing 

3.1 Overrelaxed simulated annealing algorithm 

As has been pointed out above the abelian projection procedure requires (partial) gauge 
fixing. The differential gauge condition Eq. (|3|), supplemented by the constraint of Fadeev- 
Popov operator positivity, is equivalent to finding a maximum of the lattice functional, 
Eq. (fj). Besides absolute maxima — which can in principle be degenerate even beyond 
trivial gauge transformations, such as constant ones and transformations within the max- 
imally abelian subgroup — the functional F(U) can have any number of local maxima. 



This feature resembles the Gribov problem of continuum gauge theories JTT| . While degen- 
erate absolute maxima, at least for Landau gauge, may be safely ignored ||22|| , imperfect 
gauging (where the system is stuck on some local maximum) may lead to fake physics 
results p3| . The aim of any reliable analysis must therefore be to drive the system by 



appropriate gauge procedures into local maxima as close as possible to the absolute ones. 
This would help to reduce systematic uncertainties due to gauge fixing ambiguities. A 
similar approach of gauge fixing in numerical simulations, motivated by global gauge fixing 
conditions |24|, p5| , has been advocated for Landau and Coulomb gauges in Ref. f26[ . 



Traditionally, the relaxation plus overrelaxation (RO) algorithm has been employed for 



MA gauge fixing. In Ref. |27| , we reported on the implementation of simulated annealing 
(SA), a technique that has been proven to be very useful in handling various optimization 
problems. In this SA algorithm [EM EO] the functional F(U) is regarded as a "spin action", 



S(s) = F(U°) = ±- £ Tr (s n U n ,,s n+ ,U^) (35) 

" * n,fj, 

where s n = 5^cr 3 g n resemble spin variables. The lattice fields U n ^ play the role of (almost) 
random local couplings. Maximizing the functional F{U 9 ) is equivalent to decreasing the 
auxiliary temperature T of the statistical system with partition function 

z = E ex P (^00) ■ (36) 

{»»} Vi 7 

One starts with equilibrating this spin glass at high temperature. Subsequently, T is 
decreased adiabatically. It is evident that in the limit T —>■ the system approaches its 
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ground state, i.e. the maximal value of S. The merits of SA can be phrased in the language 
of solid state physics: standard relaxation that corresponds to fast cooling might cause 
defects which are avoided in the adiabatic cooling procedure of SA. In order to enhance the 
adiabatic movement of the spin variables through phase space we complement simulated 
annealing with overrelaxation (OSA). 
Our procedure consists of three steps: 

1. Thermalization at T = 2.5. 

2. Gradual decreasing of T down to T = 0.01. 

3. Final maximization by means of the RO algorithm. 

In steps 1 and 2 an overrelaxation transformation is performed at six consecutive lattice 
sites and heatbath is applied to the seventh. Within step 2, every time when the heatbath 
update is applied to a site, the temperature is lowered by a quantum ST. For our large 
volume studies (V = 32 4 ), a variant of this algorithm (suitable for parallel computers like 
the Connection Machine CM-5 where sites are visited in lexicographical ordering within 
subcubes of 2 4 sites each) has been employed. The combined effect of local overrelaxation 
plus local temperature reduction is to cut the number of cooling sweeps while remaining 
close to equilibrium. 

For the initial thermalization at T = 2.5, 20 sweeps have been performed. Within 
the temperature range 2.5 > T > 0.1, 5T(T) has been tuned such that the spin action 
increased about linearly with the number of iteration sweeps. This has been realized by 
subdividing this range into 24 intervals of width AT = 0.1. 

The corresponding differences of the action AS(T) = S(T) — S(T — AT) have been 
computed on equilibrated configurations and were found to be very stable against statisti- 
cal fluctuations among different Monte Carlo (MC) configurations. We found practically 
no volume dependence and only a moderate impact from variations of the gauge cou- 
pling, (3. The number of sweeps (out of a fixed total number) to be performed within 
each interval (T — AT, T] was chosen to be proportional to AS(T) and, subsequently, the 
corresponding value of 5T(T) has been determined. Within the region 0.1 > T > 0.01, 
50 additional sweeps have been performed. Finally, the RO algorithm has been applied 
till a convergence criterion was satisfied^. 



2 The iterations have been stopped as soon as all rotations of s n among the lattice sites within a 
maximization sweep were equal to identity within single precision numerical accuracy. 
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3.2 Tuning the gauge fixing algorithm 

Prior to the large volume production runs, we compared the RO gauge fixing algorithm 
with three different OSA variants, using a 12 4 lattice at (3 — 2.43. The OSA variants 
employ three different cooling schedules at step 2 of our procedure: 250 (OSAl), 500 
(OSA2) and 1000 (OSA3) total sweeps. The standard procedure (RO) has been applied 
with identical convergence criterion. We collected 30 statistically independent equilibrated 
MC configurations and produced 10 random gauge copies from each of them as inputs for 
the four algorithms. 

Consider a gauge dependent abelian quantity A. In the following, we denote A to 
be the average over gauge copies and (A) to be the statistical average. In Table [l|, 
we compare the output of the four algorithms for various quantities. 5 2 d = (F 2 — F ) 
denotes the scatter of the maximized value of the functional among gauge copies. An ideal 
algorithm would always yield F = F max , i.e. 6 sd = 0. p mon = 1/(4V) k n>IJ ) denotes 
the so-called monopole density. Wy are abelian i x j Wilson loops and K ab is the abelian 
string tension in lattice units. 

The comparison of (F) and 5 s d reveals that anyone of our OSA schedules is superior 
to the standard RO algorithm, with the longest schedule yielding the best results. The 
OSA2 and OSA3 algorithms perform equally on the level of statistical errors. In terms 
of total computer time spent with a scalar code, OSA2 is about a factor two slower than 
RO. OSAl, still being an improvement, is only 30 % slower than RO. In our parallel 
implementation on a Connection Machine CM-5, OSAl was found to run even slightly 
faster than RO. 

The lesson from Table |] is that physical results can be substantially corrupted by 
inadequate gauge fixing. The functional F turns out to be correlated with Wilson loops 
and anti-correlated with the monopole number and abelian string tension. 

By applying state-of-the-art smearing techniques on the spatial transporters of abelian 
Wilson loops, we have been able to compute the abelian ground state potential. In 
Fig. m these potentials are displayed for the RO (used by all previous authors for MA 
gauge fixing) and the OSA3 algorithms with only one gauge copy used in both cases. 
In addition, the copy with largest functional among all OSA3 copies has been chosen 
on each configuration as our best estimate of the "true" maximum (circles). From this 
potential we obtain K^ est = .0478(38) as an estimate of the abelian string tension^]. From 
Fig. [I] it is evident that this value will be overestimated by about 30% by use of the 
standard procedure (see also Table [[]). Even our most expensive algorithm OSA3 yields 
a value that is off by about one statistical standard deviation from our best estimate 
(i^ SA3 = .0536(30)). 

3 Within this section, all errors have been obtained by the jackknife procedure. 
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In order to study the scaling properties, the investigations have been partly repeated on 
a lattice with nearly equal physical size but smaller lattice resolution (16 4 at (5 — 2.5115). 
The OSA2 and RO algorithms have been applied to 10 copies on 20 configurations, each. 
The qualitative properties are seen to be the same. Again, the string tension obtained 
after application of the RO algorithm is drastically overestimated: K^ Q = .0362(15) 
versus K^ est = .0305(19). Since the non-abelian string tension at this (3 value turned out 
to be K = .0325(12), this difference is relevant to the physical conclusions drawn. 

We conclude that the quality of MA gauge fixing (in terms of the value of F(U 9 ) 
and the scatter of results among different gauge copies) can be significantly improved by 
applying the OSA algorithm, without any cost in computer time. Our systematic study 
of non gauge invariant quantities revealed that such improvement is in fact mandatory: 
for a reliable extraction of the abelian potential and other observables one must carefully 
eliminate biases from incomplete gauge fixing. An extension of the method to other gauges 
(e.g. Landau gauge) is straight forward. 

4 Biases from Gribov copies 
4.1 Simulation technicalities 

Our main simulations have been performed on 32 4 lattices at /3 — 2.5115. As a first 
step, test runs with different OSA cooling schedules (as explained in the previous section) 
have been performed on four gauge copies, generated from two thermalized configurations 
with the above simulation parameters. We have chosen the following schedules for the 
temperature degrading sweeps: N s = 250 + 50, N s = 500 + 50, N s = 1000 + 100, iV s = 
2000 + 100, and N s = 5000 + 200. The first numbers denote the sweeps spent within 
the linear region (0.1 < T < 2.5). The latter numbers correspond to the sweeps applied 
within the interval 0.01 < T < 0.1. Finally, direct maximization sweeps (T = 0) have 
been applied. Depending on the configuration and cooling schedule typically 10-100 such 
steps had to be performed until the rest vector criterion was satisfied. 

The hysteresis curves of the spin action, F, as a function of the temperature, T, are 
displayed in Fig. |[ The upmost curve corresponds to the longest schedule. The scatter 
within each curve (which is not visible on the scale of the figure) is due to the four gauge 
copies generated and indicative for the (small) statistical uncertainty. The differences 
between the curves indicate that even for N s > 2000 thermal equilibrium is not yet 
reached. 

For our final run, we decided to apply 1100 gauge fixing sweeps as a compromise 
between effort and outcome. This choice allows to generate a few local maxima of F(U) 
on each configuration and subsequently to select the best out of them. As will be described 
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below, the availability of several local maxima enables us to estimate systematic biases 
due to incomplete gauge fixing. 

4.2 Error estimates 

In order to estimate systematic errors induced by an incomplete gauge fixing procedure, 
we have generated N = 20 random gauge copies on each of N c = 30 gauge configurations. 
Subsequently, these copies have been fixed to the MA gauge. Each abelian configuration 
(i.e. each gauge copy) Cj(i) = {0^' l '{n) : n — 1, . . . , 4, n G V} is labeled in the following 
way: j runs from one to the number of gauge copies, N, while i runs from 1 to the 
number of Monte Carlo generated SU{2) gauge configurations, N c . On each of these 
copies, abelian quantities A(Cj(i)) are measured, where A denotes either a (smeared) 
Wilson loop, the plaquette or the monopole number. In addition, the values of the gauge 
fixing functional F(Cj(i)) are stored. All gauge copies (on a given configuration) are 
sorted by the value of this functional: 

F(d) < F(C 2 ) < ■ ■ ■ < F(C N ) . (37) 

We are now prepared to investigate, which value of A we would have obtained on the 
"best" out of m < N copies. To answer this question one has to select m random copies 
out of the iV copies that have been generated in total and subsequently extract A from 
the copy with largest functional F, i.e. largest index i. Averaging over all possible choices 
yields 

A(Cj) (38) 

as the "average best copy" expectation of A on subsets of size m where the gauge con- 
figuration index has been omitted. As expected, the above formula corresponds to the 
average over all copies for the special case m = 1 while for m = N one obtains the value 
of A on the "best" copy. Quantities like the abelian potential can be computed from 
averages of these A m 's over the N c configurations. The situation is visualized in Fig. [3] 
for the gauge functional F itself. Throughout the rest of this paper, all statistical errors 
have been computed by the bootstrap procedure. 

We wish to determine the expected deviation of the estimate (A m ) = Ara(0 
from the true value (A^), related to an absolute maximum. To estimate this bias we 
make the following assumptions: 

• The absolute maximumF] is unique. 

The term "absolute maximum" is to be understood modulo trivial degeneracies, due to constant 



Am \ m ) \£w-™) 
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• The algorithm is in principle able to reach the absolute maximum F g iobai> i- e - Foo = 

-^global • 

• In the neighborhood of the absolute maximum, (A m ) approaches (A^) as a mono- 
tonous function of (F^ — F m ). 

The first assumption is supported by Zwanziger's proof of the non-degeneracy of absolute 
maxima within the interior of the fundamental modular region in Landau gauge on the 



lattice while the third assumption is supported by numerical evidence as all our 
observables exhibit strong correlations with the average value of the functional. 

The difference, A^ = {A^ — A m ), is to be seen as the very bias on A m from incomplete 
gauge fixing. The statistical uncertainty on this bias, on the other hand, is nothing else 
but the systematic error on our final result. As to its m-dependence, we start from the 
ansatz 

A^ = d exp(-dim) + • ■ ■ (39) 

for the (large m) asymptotic behavior, to be tested against the data. Accordingly, we can 
fit our data to the form 



(A m ) = (Aoo) - ci exp(-rf 1 m) + • • • (40) 

where (A^), C\ and d\ are free parameters. The statistical errors on (A m ) imply a large 
uncertainty on {Aoo). However, due to strong correlations among the data, the differences 
A^ can be obtained quite accurately. In view of our limited statistics for the study of 
gauge fixing ambiguities (N c = 30) we have not applied full correlated fits. Nonetheless, 
some of the correlations have been taken into account by fitting the differences 

A A (N, m) = {A N - A m ) = Ai ( e C*-»0* - l) , A^ = c x e- N * (41) 

rather than following Eq. ^0|. The result of such a fit (from m = 5 onwards) is visualized 
in Fig. |J On the basis of this fit the bias on A m can be traced back into the region of 
small m, in the form A^ = A^ + A A (iV, m). 

Results from the above procedure for the functional, F, the abelian plaquette action, 
5 = 1 — {W ab (l, 1)), and the monopole density p mon for m = 1, 10, 20 as well as for the 
extrapolated value (m = oo) are compiled in Table ^. We find that the proposed error 
analysis is a powerful tool to obtain reliable estimates for biases. For instance, when 

gauge transformations, gauge transformations within the unfixed diagonal U(l) subgroup and the 2 4 
degenerate maxima from Zi center group transformations within hyperplanes perpendicular to the four 
possible lattice orientations, which only affect expectation values of Polyakov line-like operators but leave 
the spectrum invariant. 
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selecting the best out of 10 copies generated by our OSA algorithm, we find AFi = 
0.000055(10), AS W = -0.000037(30), and Ap morijl0 = -0.000015(15). 

We would like to emphasize that, by computing the biases, we have found a way to 
extrapolate values, obtained on local maxima, to an absolute maximum. The accuracy of 
all computations is limited by the statistical error on the biases. 

4.3 Application to the abelian potential 

The abelian potential, V^f (R), as well as the abelian string tension, K^, have been com- 
puted for various m < N, by use of the method described in Section 5.1. We find 
the abelian string tension, to be anticorrelated with the gauge fixing functional 
(K^ < i^^ ) _ 1 ). The potential values themselves exhibit a systematic, but statistically 
insignificant, drift. 

The correlation between the fitted abelian string tension and m can be read off 
from Fig. [|. In analogy to the extrapolation method discussed in the previous sec- 
tion, A K (N, m) is fitted to an exponential ansatz. As a result, the systematic bias, 
A^ = A^r — A K (N, m) is obtained as a function of m (Fig. ^]). The statistical errors on 
N c = 30 configurations^] (solid line) are indicated in the figure, as well as the error on 
the final statistical ensemble of 108 configurations, obtained on the best out of 10 gauge 
copies (horizontal dashed line). 

For the final production runs we decided to choose m = 10 as a reasonable compromise. 
For this m-value the expected statistical error on the string tension matches both, the 
size and the uncertainty of its bias. 

Contrary to the case of the string tension (Fig. ^j, we find no statistically significant 
bias on the values of the Coulomb coefficient e and the self energy Vq which are dominated 
by the short range part of the potential. 

Results on the abelian string tension from different numbers of gauge copies are col- 
lected in Table Q (last row). The bias on 10 gauge copies is A^ = —0.00078(72) while 
the expected bias on one gauge copy would have been Af = —0.0021(9). 

5 Physics results 

Our main measurements have been performed on the same lattice volume and (3 value 
as the investigation of systematic gauge fixing errors, presented in the previous section 
(V = 32 4 , f3 = 2.5115). This enables us (a) to correct the results for the estimated 

5 The statistical errors come out to decrease with decreasing m due to the procedure of combinatoric 
averaging, Eq. [38| 
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biases and (b) to include systematic uncertainties into the final errors. The spatial lattice 
extent corresponds to 2.7 fm in physical units where the scale has been obtained from the 
value k = Ka~ 2 = (440 MeV) 2 for the string tension. This ensures finite size effects on 
the potentials to be negligible. The (3 value was chosen sufficiently large to be within the 
scaling region, and in respect to future finite temperature studies (The critical temperature 
comes out to be T c = (8a) _1 |j30| ). 



The hybrid-overrelaxed algorithm [31] with Fabricius-Haan heatbath sweeps |32|] has 
been applied to update the gauge fields. Subsequent configurations are separated by 200 
such sweeps and have been found to be statistically independent. Abelian projection 
and measurements have been performed on 108 such gauge configurations. The follow- 
ing abelian observables (i.e. quantities expressed in terms of abelian gauge fields, n ^) 
have been investigated: charge one and two potentials and string tensions, photon and 
monopole contributions to the abelian string tension and the abelian monopole density. 



5.1 Abelian and non-abelian static potentials 



We have computed the abelian potential on 108 configurations as well as the non-abelian 
potential on 644 configurations by use of one and the same analysis method to allow for 
direct comparison of results. 

Our final results on the abelian potential have been obtained on the best out of 10 
gauge copies. The systematic biases are estimated from 30 configurations with 20 gauge 
copies on each. Both, abelian and non-abelian potentials have been obtained from Wilson 
loops with smeared spatial paths in order to enhance the overlap of the QQ creation 



operator with the QQ ground state | 19fl . To reduce statistical fluctuations, we have 
analytically integrated out temporal links of the non-abelian Wilson loops p3| . 
We found good results by iteratively applying the smearing procedure, 



On, 



arg 



aexp 



tie. 



n,j. 



+ ex p (* {p 

k+j 



ii, k 



+ 9 



n+k,j 



n+j,k 



(42) 



to spatial link angles, i.e. by substituting the corresponding £7(1) element by a linear 
combination of the previous one and the sum of the four spatial staples, enclosing it. We 
have chosen the parameter value a = 1 and 150 iterations. Rectangular Wilson loops, 
constructed from such smeared spatial links, can be decomposed into a linear combination 
of various loops with fixed corners but different spatial connections. For extraction of the 
non-abelian potential, we applied a similar procedure on the SU(2) link variables with 
the parameter value a = 2. 
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For large temporal extent, T, the potential can be extracted from the asymptotic 
expectation, 

(W(R, T)) = C(R) exp (— V(R)T) + • • • , (43) 

where C(R) denote the ground state overlaps. For finite values of T, we define the follow- 
ing approximants to overlaps and potentials, which will monotonously decrease towards 
their asymptotic values, 

V T (R) = log L^^^L\ , C T (K) = (W(R,T))exp(V T (R)T) . (44) 

By comparing results on the potential approximants (abelian and non-abelian) be- 
tween the on-axis direction with those, obtained on five different off-axis directions, we 
find 5*0(3) rotational invariance to be restored (within statistical accuracy) for R > 3. In 
both cases, the data for R > 2^/3 are well described (for T > 3) by the parametrization, 

V(R) = V + KR-^ . (45) 

By fitting different T-approximants to the potential with this parametrization, we obtain 
approximants to the string tension, Kt- By demanding plateaus of Ct{R), Vt{R) and Kt 
for T > T m i n , we find T m j n = 4 and T m j n = 6 for the abelian and non-abelian potentials, 
respectively. The different onset of asymptotics is due to superior overlaps with the 
ground state in case of the smeared abelian Wilson loops. The self energy term, VQ b 
is much smaller than its abelian counterpart Vq, resulting in larger numerical values of 
the corresponding abelian Wilson loops. These two effects are among the reasons for 
reduced statistical errors on the abelian potential and fit parameters. We take the T m i n 
approximants to potential values and fit parameters as our asymptotic results. To avoid 
systematic effects from the fit range creeping into the comparative interpretation of results, 
we select a universal i?-range for all fits, 2y/3 < R < 16. This provides us with 45 on- 
and off-axis data points. 

Our results on the fit parameters are in qualitative agreement with previous publica- 
tions [§]. In particular, the self energy and Coulomb coefficient e ab come out to be 
by more than a factor two smaller than their non-abelian counterparts while the abelian 
string tension is found to be close to the non-abelian one (see Table | and Fig. |7|). As 
pointed out above, all systematic uncertainties are understood and under control in the 
present investigation, in particular the approach to the T — > oo limit. The biases due 
to gauge fixing ambiguities have been neglected in previous studies but turn out to be 
important as demonstrated in Sections 3 and 4. 

The value of the abelian string tension comes out to be K ab = 0.0305(3) on the best out 
of 10 OSA gauge copies, which nicely agrees with the value K ab = 0.0305(19), as obtained 
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on a 16 4 lattice (Section 3) at the same (3 value. By including bias and uncertainty of the 
bias, we end up with K ab = 0.0297(8), where the error includes the systematic uncertainty. 
This amounts to the ratio K ab /K = 0.92(4) (see Table g). 

Given the precision of our data and analysis tools, we find a significant deviation 
between the string tension of SU(2) gauge theory and its MA content at (5 — 2.5115. 
Further measurements at different lattice spacings are required to decide whether this ratio 
approaches unity in the continuum limit, as expected from abelian dominance. However, 
for the time being, at finite lattice spacing, the string tension from the abelian projected 
theory appears to be definitely smaller than its non-abelian counterpart. 

5.2 Decomposition of the abelian potential 

We have tried two different approaches to disentangle monopole and photon contributions 
to the potential. 

The first method rests on the determination of the monopole part from elementary 
monopole currents, k n ^. The potential estimators V^ on (i?), extracted from (W mon (R, T)) 
(see Eq. (0)) at fixed R have been observed to increase with T, which means that the 
coefficients C'(R) of the decomposition, 

(W mon {R, T)) = C{R) exp {-V rnon {R)T) + C'{R) exp {-V' mon {R)T) + ■■■ , (46) 

are not necessarily positive. This unpleasant feature, in conjunction with the requirement 
of T 3> R makes this procedure of extracting V mon (R) unreliable. 

As a way out, we started from an alternative representation of the factorization prop- 
erty to W ab (Eqs. fl2"6| ) and (P7|)), which allows decoupling of excited states through 
smearing of spatial links 9™° n . In this second approach , the abelian configuration has 
been fixed for technical reasons to Landau gauge, 

<9 M sin(^)=0 , (47) 

prior to evaluation of 6*™° n (Eq. (|26|) ). This reduces the number of plaquettes with nonzero 
m n ^ v and, thus, computer time by more than one order of magnitude. It should be noted 
that — though transforms under £7(1) gauge transformations Eq. (]5|) — \y mon J v as 
well as W ph are gauge invariant quantities. 

After evaluation of 6'™°™ and 9' = #„ iM — smearing has been applied and the 

corresponding potentials have been extracted. For small R, the results as obtained from 
the two methods are found to agree, while at large R, plateaus in T could only be es- 
tablished for potential estimates extracted from smeared monopole and photon Wilson 
loops. 
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As pointed out in Section 2, the full factorization ansatz may include finite volume 
contributions. However, no statistically significant such effects have been found by com- 
parison of unsmeared W mon ^ v with W mon (calculated by use of elementary monopole 
currents k n ^ via Eq. ([18]) for all realized R and T values, on a configuration by config- 
uration basis). Thus, we conclude that finite size effects can indeed be neglected on our 
lattice volume. 

The resulting potentials as well as V ab are displayed in Fig. || and the parameter values 
as obtained from fits according to Eq. (flop are quoted in Table f|. The corresponding values 
for the abelian potential are included as well, where we have omitted the bias from gauge 
fixing ambiguities, to allow for a direct comparison. 

One can see from Fig. [| that the photon part V ph does not contribute to the string 
tension. Therefore, in MA projection only the monopole part of the abelian gauge fields, 
9™™, gives rise to the flux tube. A comparison of V ab with V mon ^ v + V ph reveals 
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qualitative agreement, implying the approximate validity of the factorization ansatz. The 
string tension of the monopole contribution amounts to (95±1)% of the full abelian string 
tension. The approximate decomposition of the abelian static potential into monopole 
and photon parts gives evidence for the interaction term between monopoles and photons 
within the corresponding effective action to be weak. 

We have also attempted to fit the photon part of the potential to the ansatz, 

4-7T 

V ph (R) = V scU -fG L=32 (R) , G L (R) = — £e* kR D(k,0) (48) 

v k^O 

for various fit ranges. Gx(R) is the Coulomb potential on the lattice which approaches 
1/R in the infinite volume limit (L — > oo) for large (lattice) R. Lattice artefacts turned out 
to be well parameterized by this functional form. As a result we quote V^ e if = 0.2513(3) 
and / = 0.130(30), obtained on the fit range 3 < R < 12.12. Data and fit curve are 
visualized in Fig. |9|. The fit range turned out to be correlated with the parameter value 
/ in so far as / tended to be larger if large R values had been included and smaller 
for small R values. This is also evident from the figure and might be interpreted as a 
relict of asymptotic freedom within the abelian projected gauge theory. From tree level 
perturbation theory, one might expect / = V^ e if/Gi = 3 2 (0) ~ 0.0814, which is smaller than 
the fitted / quoted above, in accordance with a running coupling interpretation. 

The monopole contribution to the abelian potential V{^ 2 has been extracted from 
extended monopoles of size I = 2 as well, by use of the method, introduced in Section 2 
(Eq. (|24])). By computing the potential from the corresponding Wilson loops (which has 



been done in Ref. \W\ the first time), we found the approximants to decrease monotonously 
in T. The value of the monopole string tension turned out to be slightly smaller than 
the one extracted from elementary monopoles (KJ^ 1 = 0.0271(3) instead of K mon = 
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0.0290(3), c.f. Table |4j), which might indicate that a small portion of the string tension is 
due to monopole structures of smaller extent. 

The new method of computing the monopole contribution to the potential by com- 
bining Eqs. (|26|) and (27]) with the smearing method has the advantage of large ground 
state overlaps. Also, the potential can easily be computed for off-axis points while the 
blocking method even reduces the number of on-axis points at which the potential can be 
measured by a factor two. Notice, that V mon ^ s obtained by use of the smearing method 
is identical to V mon + V^ s , as computed in the traditional way from elementary monopole 
currents, while monopole structures on the scale of a lattice spacing are neglected in the 
blocking method. 

5.3 Charge two case 

In order to extract the charge two potential V ab,2 (R) (Section |2.4| ), smeared charge two 
Wilson loops, W ab,2 (R, T) have been evaluated. The data for the fundamental and charge 
two abelian potentials are displayed in Fig. (with the fitted self energies being sub- 
tracted). 

Flux tube models lead to the expectation K ad i / K = 8/3 for the ratio of the adjoint 
over the the fundamental string tension^. This value was qualitatively supported by 



numerical data [33, ^5|. However, recent results [p6 , [37| indicate that at intermediate to 
large distances the above ratio tends to be somewhat smaller than 8/3. We suggest the 
slope of the charge two potential, K ab ' 2 to constitute the abelian projection counterpart 
to the adjoint string tension and thus expect the ratio K ab ' 2 / K to agree with results on 
K i I K. However, the lack of high precision data on the latter ratio prevents us from a 
quantitative test of this assumption. From our data, we find the value K ab,2 /K ab = 2.23(5) 
significantly smaller than 8/3. 

For small R, perturbation theory yields the same 8/3 ratio between the adjoint and 
the fundamental potential, which differs from the expected ratio between abelian charge 
two and charge one potentials, where one naively would expect a factor 2 2 = 4. We 
indeed obtain V^ 6 ' 2 /Vg* = 3.96(7) and e ab ' 2 /e ab = 4.0(5) for the fit parameters that are 
sensitive to short range physics. We have included the curve 8/3K ab — 4e ab /R into Fig. pi], 
which corresponds to SU(2) Casimir scaling for the linear part of the potential and to 
the perturbative U(l) expectation for its Coulomb part. 



6 Of course, the adjoint string tension is only an effective quantity since, at large distance, string 
breaking is expected to set, caused by screening from the creation of glueball pairs. 
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5.4 Monopole density 



Monopole densities for elementary as well as for extended monopole currents of size I = 2 
have been evaluated. The density of monopoles of size I in physical units is defined as[| 



n (l),ph 

r mon 



AV ■ {la} 3 



E(lCl) 



(49) 



n,)i 



We obtain the following results for the monopole densities, converted into units of the 
string tension as measured in the present investigation: 



r mon 
#3/2 



1.962(4) 



r mon 



1.269(2) 



(50) 



The value for elementary monopole currents has been corrected by its systematic bias 
from incomplete gauge fixing. 

In order to relate our results to those obtained in previous publications, we compare 
our value on the density of elementary monopoles with the value p^ on / K 3 ^ 2 = 2.11(2) from 
Ref. [Q for a 16 4 lattice at /3 — 2.5, which is close to our coupling constant, f3 = 2.5115. 
Note, that we have rescaled the result of Ref. |38] into units of the string tension as 
ill . The result of Ref. turns out to be consistent within errors 
I7L H0[ while our value is significantly smaller. Since 



obtained in Ref. 



with results from other authors 39 



finite size effects are negligible for the lattice extents under consideration, the difference 
seems to be due to our improved OSA gauge fixing algorithm. A similar discrepancy is 
observed for extended monopoles. In this case our value, Eq. (|50|) , should be compared 
with p^l h /K 3 l 2 = 1.32, as obtained on a 24 4 lattice at /3 = 2.5 by means of the standard 



OR gauge fixing procedure [41 



6 Summary and conclusions 

Let us summarize our main results and conclude: 

• The present study is based on a self contained and self consistent analysis on the 
largest lattice volume that has been studied so far for this kind of simulations with 
comparatively high statistics. A systematic error analysis has been carried out for 
the first time in this context. 

7 The factor I 3 appears in the denominator of Eq. (pL9|), due to averaging over the I 3 possible blocked 
sublattices. 
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• To obtain reliable results in the abelian projected theory with MA gauge condition, 
one has to investigate and control the uncertainty that is inevitably introduced by 
the incomplete gauge fixing of numerical practice. The OSA algorithm has been 
shown to be a powerful tool for gauge fixing. A method for estimation of residual 
uncertainties is proposed. 

• We have found K ab to be (8 ± 4)% smaller than the non- abelian string tension at 
P = 2.5115. 

• Our investigation of the decomposition of the static abelian potential into monopole 
and photon parts confirms earlier observations at a higher confidence level. By 
applying a new method for extracting the monopole contribution to the potential, 
we have been able to extract the corresponding potentials and fit parameters reliably. 
The factorization has been found to work qualitatively, the monopole contribution 
accounting for the string tension within a margin of 5%. 

• We have calculated the abelian projection approximation for the adjoint string ten- 
sion. Our result for ratio K ab,2 /K ab is in qualitative agreement with numerical data 
on its non- abelian counterpart. 

We believe that further computations at different (3 values, provided all sources of 
errors are kept under control as in the present paper, will answer the question whether 
the abelian projected theory exactly reproduces the large distance behavior of the full 
theory in the continuum limit. 
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Note added in proof 



After completion of the present paper, a preprint by G. Poulis [42] has been received, in 
which the author has demonstrated that one should expect K ab > 2 and K ad i to coincide to 
the same extent as K ab and K agree by use of reasonable approximations. 
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Table 1: Comparison of gauge fixing algorithms. 



RO OSA1 OSA2 OSA3 



(F) 


0.7370(2) 
18 • 10~ 5 


0.7383(2) 
13 • 10~ 5 


0.7387(2) 


0.7390(2) 
8 • 10~ 5 




10 • 10~ 5 


Pmon 


0.0218(2) 


0.0209(3) 


0.0207(3) 


0.0204(3) 


W n 


0.7702(5) 


0.7716(5) 


0.7720(5) 


0.7723(5) 




0.085(1) 


0.090(1) 


0.091(1) 


0.092(2) 


R ab 


0.063(3) 


0.057(3) 


0.055(3) 


0.054(3) 



Table 2: Dependence of results on the number of gauge copies, m, at f3 — 2.5115, V = 32 4 , 
N c = 30. 





171=1 


m — 10 


m = 20 


m = oo 


(F) 


0.752641(34) 


0.752789(35) 


0.752822(36) 


0.752845(37) 


S 


0.202508(60) 


0.202299(64) 


0.202267(64) 


0.202261(66) 


Pmon 


0.011639(26) 


0.011528(30) 


0.011515(32) 


0.011513(33) 


K ab 


0.0325(11) 


0.0311(13) 


0.0306(16) 


0.0303(17) 



Table 3: Fit parameters for the static potentials in abelian projected SU(2) (corrected 
by the estimated bias due to incomplete gauge fixing) and the potential of full 577(2). 
Errors are systematic and statistical. 





AP SU{2) 


SU(2) 


e 


0.095(11) 


0.252(24) 


V 


0.240(4) 


0.545(10) 


K 


0.0297(8) 


0.0325(12) 



Table 4: Test of factorization of the abelian potential into monopole and photon parts. 
Biases from incomplete gauge fixing have been omitted throughout the table. 





yab 


ymonjv _|_ yph 


ymonjv 


l/mon 
v l=2 


yph 


e 


0.095(6) 


0.068(5) 


-0.056(5) 


0.019(4) 


0.124(2) 


V 


0.240(3) 


0.232(2) 


-0.029(2) 


0.010(2) 


0.261(1) 


K 


0.0305(3) 


0.0291(3) 


0.0290(3) 


0.0271(3) 


0.00007(4) 
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Table 5: Static potential parameters for the abelian projections of the fundamental (q = 1) 



and adjoint (q = 2) representations of SU(2). 





AP SU{2) (q = 1) 


AP SU{2) (q = 2) 


(q = 2)/{q = 1) 


e 


0.095(6) 


0.376(36) 


3.98(48) 


V 


0.240(2) 


0.950(13) 


3.96(7) 


K 


0.0305(3) 


0.0682(11) 


2.23(5) 



0.6 



0.5 



£ 0.4 - 



0.3 - 



i i 

RO — A — ' 
OSA3 — a — ' 
best e i 



0.2 - * 



16 



[] 

O 



1 2 3 4 5 6 

R 



Figure 1: Abelian potentials from RO (triangles), OSA3 (squares) and "best" copy (cir- 
cles). 
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Figure 2: The spin functional F as a function of the temperature T for various cooling 
schedules. The curves correspond to N s = 5200, N s = 2100, N s = 1100, iV s = 550, and 
N s = 300, respectively. 
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Figure 3: The maximized value of the spin functional F as a function of the number of 
local maxima generated on each configuration. The solid line (with error band) represents 
our extrapolated value (m — > oo). 
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Figure 4: The differences, A F (N,m) = F N — F m and an exponential fit. 
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Figure 5: The abelian string tension versus the number of MA gauge copies m. The solid 
line (with error band) represents our extrapolated value (m — > oo). 



30 




Figure 6: Differences between the extrapolated value (from the exponential fit) and 
K ab , obtained on a finite number of gauge copies, m. The solid line denotes the statistical 
uncertainty on K ab from 30 gauge configurations. The horizontal dashed line is the 
statistical uncertainty on 108 configurations with m = 10. 
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Figure 8: The abelian potential (diamonds) in comparison with the photon contribution 
(squares), the monopole (plus finite volume) contribution (crosses) and the sum of these 
two parts (triangles). Notice, that no (self energy) constants have been subtracted from 
any of the data sets. 
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Figure 9: The photon contribution to the abelian potential. The dashed vertical lines 
indicate the fit range (solid curve and squares). 
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Figure 10: The static abelian charge two potential V ab,2 (R) in comparison to V ab (R). 
In addition, the expectation V ab > 2 (R) - V^' 2 « 8/3K ab R - Ae ab /R is included (upmost 
dashed line). 
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